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ABSTRACT 

It is known from numerical simulations that thermal conduction along magnetic field lines 
plays an important role in the evolution of the kink instability in coronal loops. This study 
presents the observational signatures of the kink instability in long coronal loops when parallel 
thermal conduction is included. The 3D nonlinear magnetohydrodynamic equations are solved 
numerically to simulate the evolution of a coronal loop that is initially in an unstable equilibrium. 
The loop has length 80 Mm, width 8 Mm and an initial maximum twist of $ = 11.57T, where 
$ is a function of the radius. The initial loop parameters are obtained from a highly twisted 
loop observed in the TRACE 171 A waveband. Synthetic observables are generated from the 
data. These observables include spatial and temporal averaging to account for the resolution 
and exposure times of TRACE images. Parallel thermal conduction reduces the maximum local 
temperature by up to an order of magnitude. This means that different spectral lines are formed 
and different internal loop structures are visible with or without the inclusion of thermal conduc- 
tion. However, the response functions sample a broad range of temperatures. The result is that 
the inclusion of parallel thermal conductivity does not have as large an impact on observational 
signatures as the order of magnitude reduction in the maximum temperature would suggest; the 
net effect is a blurring of internal features of the loop structure. 

Subject headings: instabilities — magnetic fields — magnetohydrodynamics (MHD) — Sun: corona 



1. Introduction 

The magnetohydrodynamic (MHD) kink in- 
stability in cylindrical geometry serves as a 
first approximation for this instability in coro- 
nal loops. Its evolution has been studied nu- 



merically for more than a decade dMikic et al 



1990HBatv fc Hevvaertslll996l ; lLionello et al.lll 998; 



Arber et al.1 1999). Previous work has studied 
the onset of the kink instability by rotating a 
flux tube's footpoints un til the twist increased 
beyond the critical value ([Galsgaard fc Nordlund 
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1997[ iGerrard et al.ll2002l ). Curvature effects have 
also been incorporated into the model by study- 
ing a kink-unstable flux tube, curved in a half- 
torus with, its footpoints anchored in the same 
plane dTorok et all 120041 iTorok fc Kliernl 120051: 
Gerrard et al.ll2004l) . 

In another series of studies, the ideal MHD kink 
instability has been evoked as a t rigger for recon- 



necti o n occurring in coronal lo ops (Browning & Va n der 
20031: iBrowning et al.l 120081: Hood et alj l2009t 



2.1. Physical model 

Nonlinear three-dimensional simulations are 
performed using the MHD Lagrangian-rema p 
code, Lare3d, as described bv lArber et al. ( 200ll ). 
It solves the resistive MHD equations for a fully 
ionised plasma, with a heat flux included in the 
energy equation. A full description o f the physical 
mode l is given in iBotha et al.1 ([201 II) . 

i nermal conduction is included along the mag- 



Bareford et al.l I2010I ). A cylinder is initialised 



with a twisted magnetic field profile that is unsta- 
ble, and is t hen perturbed to ev olve into the kink 
instability. IBotha et al. ( 2011 ) have shown that 
by including thermal conduction along magnetic 
field lines in these models, the maximum tem- 
perature obtained during reconnection events is 
lowered by an order of magnitude, activating dif- 
ferent spectral lines than when no parallel thermal 
conduction is included. 

In this paper the observational signatures of 
the MHD kink instability, with thermal con- 
duction included, are presented and compared 
with a simulation without thermal conduction. 
A kink-unstable cylinder is evolved and the nu- 
merical results are co mpared to observations of a 
twisted coronal loop ( Srivastava et al. 2010t ). In 
order to facilitate the comparison between the 
observed and forward modelled loops, the initial 
conditions are o b tained from the observations by 
Srivastava et al. (l2010h . The numerical results 
are filtered through the TRACE temperature re- 
sponse function before synthetic images in the 171 
A band are presented of spatially and temporally 
averag ed line of sight intensity measurements, sim- 
ilar to lHavnes fc Arberl (|2007t) . 

The paper is divided into two main parts. In the 
first part the model is described, which includes a 
short description of the physical model, the initial 
conditions of the numerical simulation and how 
the synthetic images were obtained from the nu- 
merical data. In the second part the numerical 
results are discussed and the synthetic intensity 
images are generated from the data. 

2. Model 

The coronal loop is initialised as a straight 
twisted cylinder in a uniform background temper- 
ature and density. 



netic fi eld in the form of the classical lSpitzer-Harm 
(|1953l ). or Braginskii, conductivity with log A = 
18.4. This corresponds to the standard ther- 
mal conductivity parallel to th e magnetic fi eld of 
«H = 10- n T 5 / 2 W m- 1 K- 1 (|Priesdl2000l ). Ra- 
diative losses do not play an important role and 
are not included in the model. For the coronal val- 
ues used in the model, the radiative cooling time is 
of order hours. During the simulations reconnec- 
tion and heating occur within short lengths along 
the loop, typically one tenth of the loop length 
( Botha et al.ll2011 ). When the conductive cooling 
time is calculated for this length, it is in the order 
of minutes. 

The code contains an artificial resistivity that is 
activated only when and where the current exceeds 
a critical value. It is of the form 



T)0, 

0. 



I J I >3c 

lil <jc, 



(1) 



where rjo is the anomalous resistivity and j c = 2 
mA is the critical current. 770 is activated as 
soon as the kink instability occurs and then it 
stays active for t he duration of the simulation 
(jBotha et alj|201lh . 



2.2. Initialisation data 

The data used to obtain the initial conditions of 
the simulated coronal loop is obt ained from obser- 



vatio ns of a highly twisted loop ([Srivastava et al 
20101) . The loop was situated in AR NOAA 10960 
and was observed on 2007 June 4 between 04:43 
UT and 04:52 UT. SoHO/MDI, Hinode/SOT G- 
band (4504 A) and Hinode/SOT Ca II H (3968 A) 
were used respectively to observe the photospheric 
and chromospheric parts of the active region and 
the associated highly twisted loop, while TRACE 
171 A was used to observe its coronal part. Fig- 
ure [1] shows a closeup of approximately half of the 
observed loop system in TRACE 171 A, 6 minutes 
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after the activation of helical twist in the loop dur- 
ing the flaring process. 

From the observations it is estimated that the 
loop length is ~80 Mm and its radius in the corona 
^4 Mm. This aspect ratio of 1: 10 is within the 
observed range for coronal loops ( Klimchukl [2 000: 
Watko fc KUmchukll200rt : lAschwanden fc Boernerl 
201 lh . The width of 8 Mm is larger than the usu- 
ally observed coronal loop width of 2-4 Mm. Given 
that active region coronal loops change between 
40% and 70% in diameter between their base in the 
upper transition region o r lower corona and their 
apex (jBrooks et al.ll2007l ). we assume that the loop 
radius at its footpoints in the upper transition re- 
gion is 2 Mm. One of the footpoints of the active 
region loop was above a positive polarity sunspot 



( Srivastava et al. 2010h and it is estimated that 



the lower bound of the magnetic field strength 
at this location was approximately 470 G. The 
minimum average strength of the magnetic field 
is estimated by assuming the homogeneous dis- 
tribution of measured magnetic fluxes above the 
selected region over the positive polarity sunspot. 
However, the strength of the magnetic field may 
change at various locations above it between its 
minimum and maximum values. From this pho- 
tospheric value the ch romospheric magnet i c field 
strength is calculated. IPetrie fc Patrikeeval (|2009t ) 
found that the photospheric field is mostly verti- 
cal while the chromospheric field has no preferred 
direction. Since we are considering a magnetic 
loop, the assumption is made that most of the 
photospheric magnetic field goes to the chromo- 
sphere. However, this w ould be an upper limit. 
Kozlova fc Somovl (|200flh looked at 36 sunspots 



and found that from a height of 100 km to 1850 
km the magnetic field strength diminishes between 
0.1 and 0.7 G km -1 , with the mean field gradi- 
ent 0.35±0.05 G km -1 . Using a magnetic field 
gradient of 0.2 G km -1 and the field strength at 
the photosphere, the magnetic field strength at a 
height of 1850 km is calculated to be 80 G. If we as- 
sume that the field strength changes little between 
the upper chromosphere and the upper transition 
region, then the flux through the loop footpoint 
(radius 2 Mm) is 10 19 Mx. With the assumption 
of flux conservation along the loop, the magnetic 
field strength at the loop apex (radius 4 Mm) is 
20 G. This lead us to initialise the loop with a 
maximum field strength of 20 G inside the loop, 



while the background outside the loop contains a 
uniform field of 15 G that is parallel to the cylin- 
drical axis. 

In the numerical simulations the coronal loop is 
initialised as a uniform cylinder in force-free equi- 
libriu m and unstable to an ideal MHD kink insta- 
bility dHood et al.ll2009h . The axial twist is given 

by 



<I> = 



LBg 

rB z 



with max($) = 11.5tt, (2) 



where L is the loop length and r is its radius. The 
axial magnetic field is given by B z and the az- 
imuthal field by Bg. Both Bg and B z are func- 
tions of r. The maximum twist is at position 
r = 1 Mm, with zero twist on the loop axis and 
at its edge. The radial profile of the twist is pre- 
sented in Figure 2(b) in iHood et all (|2009h . The 
maximum number of twist (six full turns) is cho- 
sen from the observational analysis of the highly 
twisted loop system shown in Figure [T] where 
three full turns are clearly visible. The total of 
six turns is estimated by extrapolating to the in- 
visible part of the loop by assuming a symmetric 
loop shape. Max($) exceeds the stability thresh- 
old so that the loop is kink unstable. The crit- 
ical value for twist in a cylinder was found nu- 
merically to be 4.8?r (jMikic et al.lll990l). while lin 



ear t he ory predicts it to be 2 . 5tt (IGerrard et al 



20021 ). iGalseaard fc Nordlundl il997j) found that 
the critical twist angle is a function of the loop 
diameter, field strength and magnetic resistivity 
in the model. A full description of the magnetic 
field structu r e used in the simulations is given in 
Botha et al.l (|201l[) . 

Gravity is absent in the simulations. The ini- 
tial temperature and mass density are chosen to 
be uniform and constant. The mass density has a 
value of 1.67 x 10" 12 kg m" 3 (|Young et al.ll2009h 
and the temperature 0.125 MK. This temperature 
was chosen to make the evolution of the kink in- 
stability visible for the TRACE satellite. During 
the evolution of the kink instability the temper- 
ature increases loc ally where reconnection occurs 
( Botha et al. 1201 lh . Through trial and error the 
initial temperature was chosen so that these high 
temperatures will be visible using the TRACE 
171 A temperature response function. 

The evolution of the coronal loop is studied in 
Cartesian geometry (Figure [2]). The boundaries 
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in the (x, y) plane perpendicular to the loop axis 
are at ±8 Mm and reflective. Thus, there are 4 
Mm between the loop's edge and the outer bound- 
aries. This distance proved to be adequate so 
that no feedback from the boundaries influences 
the numerical results. Along the loop axis the 
boundaries are at ±40 Mm with velocities held at 
zero and the temperature fixed at the initial back- 
ground value but allowing temperature gradients. 
Hence, a heat flux across the ends of the loops 
exists. The grid resolution in (x, y, z) is given by 
128 x 128 x 256. 

2.3. Image generation 

The images generated from the numerical simu- 
lations were obtained by usin g the temperature re - 

sponse functions of TRACE ( Schriiver et al.lll999 ) 

as we ll as those of SDO / AIA (jAschwanden fc Boerner 
20111 ). Details of the li ne contributi o ns fo r the 
TRACE channels are in lHandv et all (Il999l) an d 



for the AIA channels in lO'Dwver et ail (|2010h 



The emission is calculated at every node of the 
numerical grid and then integrated along the y 
direction perpendicular to the axis of the loop, as 
indicated in Figure [5J The line of sight integral is 
given by 



9(T)p 2 dy 



(3) 



where / is the measured intensity, g(T) the tem- 
perature response function of the respective in- 
struments and p the mass density. Figure [3] 
presents the temperature response functions of 
TRACE and AIA for 171 A. This integration pro- 
duces intensity images in the (x, z) plane, which 
are then integrated over time, with the time in- 
terval determined by the exposure time of the in- 
strument. For TRACE this is taken to be 31.9 s 
and for SDO/AIA the time interval is 2.9 s. Fi- 
nally, the time integrated image is degraded by 
spatially averaging over squares of 0.375 x 0.375 
Mm 2 to compensate for the pixel resolution of 
both TRACE and SDO/AIA. The process is illus- 
trated in Figure IH using the TRACE 171 A tem- 
perature response function as given in Figure [3] 

3. Discussion 

Two simulated data sets are presented: with 
and without thermal conduction parallel to the 
magnetic field lines. All other parameters and 



initial conditions are kept identical between the 
two runs. The physical consequences due to the 
inclusion of parallel thermal conduction are dis- 
cussed bv iBotha etal] (|201ll) . In this paper the 
line of sight intensity contours are presented as 
observed through TRACE response functions. As 
such, only physical processes that are needed to 
explain the observables are included. The density 
and temperature from both data sets are used with 
TRACE response functions to generate intensity 
contours that can be compared with observational 
data. 

After initialisation, the kink unstable loop 
evolves through a linear phase that lasts for 300 
s, as is evident from Figure [S] During the non- 
linear phase the kink instability drives magnetic 
field into a current sheet, where reconnection oc- 
curs and the temperature reaches a maximum due 
to the energy released. Subsequently, thermal 
conduction along magnetic field lines transports 
heat along the magnetic field lines away from the 
points where the plasma was heated by reconnec- 
tion. The time scale for the evolution of the kink 
instability is the same with and without thermal 
conduction (Figures [5] and [B]). In the case without 
thermal conduction the heat due to reconnection 
is not conducted along magnetic field lines and 
the temperature maximum is higher. This can be 
seen in Figure [5] where the maximum temperature 
reached without parallel thermal conduction is 9 
MK, compared to a maximum of 3 MK with ther- 
mal conduction. The physical processes during the 
kink instability, with and without the rmal conduc- 
tion, a re discussed in more detail by IBotha et al 
<|201lh . 



The emission, as observed in the TRACE 171 
A band, is presented in Figures [7] and [8] for the 
cases with and without parallel thermal conduc- 
tion. Figure [3] shows that the temperature re- 
sponse function lies within the range determined 
by the minimum and maximum temperatures from 
the simulation with thermal conduction (Figures 
[5] and [6J. The higher maximum temperature 
obtained without thermal conduction (Figure [5]) 
causes the temperature response function to sam- 
ple different aspects of the evolution of the kink 
instability. However, Figure [B] show that the aver- 
age loop temperatures with and without thermal 
conduction are comparable, making the intensity 
images in Figures [7] and [8] not radically different. 
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This is due to the fact that heating occurs only at 
small localised areas, from which heat is then con- 
ducted along magnetic field lines. Figures [7] and 
[5] show that the kink instability causes the same 
structures to form in both cases. Thermal con- 
ductivity causes the temperature to spread along 
magnetic field lines, resulting in images of which 
the features are less defined when compared with 
the images generated from data without thermal 
conductivity. 

Irrespective of the inclusion or exclusion of ther- 
mal conductivity, the emission images show that 
the footpoints of the loops increase their emission 
during the nonlinear phase of the kink instability. 
The line of sight integral © is determined by the 
temperature response function as well as the mass 
density. Figure |9] gives the temperature and mass 
density profiles along the central axis towards the 
end of the simulation for the case when thermal 
conduction is included. Note that the central axis 
experiences some heating (it has higher tempera- 
tures than the average temperatures in Figure [6]) 
but the locations with maximum temperature are 
not on it (as it has lower temperatures than the 
maxima in Figure [3]). Figure |H] shows that the en- 
hanced emission is due to a density increase at the 
footpoints - and not footpoint heating. Plasma 
flows from the middle of the numerical domain, 
where current sheets form and reconnection oc- 
curs, driven by MHD ponderomotive forces gener- 
ated during the kink instability. Footpoint bright- 
ening due to compression was also observed in 
the co ronal loop simulations of lHavnes fc Arber 
(j2007h . 

One factor that determines the onset of the non- 
linear phase of the kink instability is the initial 
twist ([2]) in the coronal loop. A smaller value 
of max(<I>) increases the duration of the linear 
phase. However, it was found that once the nonlin- 
ear phase is reached, the formation of the current 
sheet, the reconnection and the thermal aftermath 
have the same duration as long as the initial twist 
exceeds the stability threshold. 

The first image in Figure [7] is sampled at 290.0 
s and shows the twisted structure towards the end 
of the linear phase of the kink instability. In the 
second image (at 321.9 s) the nonlinear phase has 
formed a current sheet and in the third image 
(at 353.8 s) thermal conduction has transported 
heat from the reconnection site along the mag- 



netic field lines. From 321.9 s the evolution of 
the internal loop structure along its length is ob- 
served for 4 minutes until the end of the numeri- 
cal run at 577.0 s. As reconnection occurs within 
the loop, magnetic fiel d lines straighten out a long 
the length of the loop (|Havnes fc Arber|[2007l) . At 
the same time heat is transported along magnetic 
field lines away from reconnection sites, so that 
the loop coo ls and its internal structure becomes 
less defined (jBotha et al.ll2011l ). 

The differences in intensity images with and 
without thermal conduction are clear when the in- 
tensity images at 449.5 s and later times in Fig- 
ures [7] and |8] are compared. Despite the order of 
magnitude reduction in the maximum tempera- 
ture when thermal conduction is included (Fig- 
ure [5]), similar looking internal loop structures 
are visible in the intensity images with and with- 
out thermal conduction. The response functions 
capture a range of temperatures (Figure [3]) and 
because the average temperatures in the coronal 
loops are of similar values (Figure [6]), the obser- 
vational signatures are less sensitive to the inclu- 
sion of parallel thermal conduction. The main ef- 
fect due to the inclusion of thermal conduction is 
the blurring of the internal features of the coronal 
loop, which is due to the efficient conduction of 
heat along magnetic field lines. In the simulation 
with thermal conduction the maximum tempera- 
ture of the magnetic structures is closer to the av- 
erage loop temperature (Figure [5]) than the maxi- 
mum temperature obtained without thermal con- 
duction. Consequently the internal loop structure 
are not as clearly delineated in the emission in- 
tensity as when no thermal conduction is present, 
which can be seen when comparing Figures [7] and 

m 

In Figures [7] and [8] it appears that the amount of 
twist towards the end of the linear phase (at 290 
s) is approximately three turns - instead of the 
six turns dictated by max("I>) in ©. The heating 
caused by the kink instability is highly localised 
and only magnetic field lines that pass through 
the heated area are visible. The physi cal processes 
associ ated with this are discussed by iBotha et al. 
(|201 lh . As the twist is a function of the radius, 



one should not expect the line of sight integral (0 
to capture the heating that occurs at max($). 

When comparing the numerical results with 
observations, it is more realistic to compare the 
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results obtained when thermal conduction is in- 
cluded in the model. Figure Q] s hows a loop seg- 
ment of the observations from ISrivastava et al 



(l201dh at a time in its evolution when the twisted 
threads are most visible in the observed images. 
The simulation images in Figure Q] shows a similar 
fine internal structure, but it is clear that a per- 
fect visual or quantitative match is not possible. 
As an example of this, Figure [7] at time 353.8 s or 
385.7 s shows the twist to be one full turn over 40 
Mm, while the observed image shows a full turn 
to be over 20 Mm. 

In this paper the results were obtained from 
only one initial equilibrium state. It is possible 
that other initial conditions may lead to end re- 
sul ts that correspond clos er with the observations 
bv ISrivastava et aL ( 2010h . This needs further in- 
vestigation. We also ignored curvature and grav- 
ity, both of which may be important during the 
evolution of the kink instability. Assuming the 
shape of the loop to be a perfect half-circle, the 
top of the loop will be 25 Mm above the transi- 
tion region. The average and maximum temper- 
atures inside the loop during the evolution of the 
kink instability are approximately 0.3 MK and 2 
MK. These temperatures give gravitational scale 
heights of approximately 20 Mm and 120 Mm. 
Thus, density stratification is important and grav- 
ity cannot be ignored in the forward modelling of 
the kink instability for this loop. The inclusion of 
both curvature and gravity should be investigated 
further. 

The observations by ISrivastava et al. ( 2010f ) 
show that after the heating event, the loop struc- 
ture cools down within five minutes to photo- 
spheric temperatures. The duration of the simu- 
lations presented here, shows 4.5 minutes of the 
nonlinear evolution of the kink instability. Dur- 
ing this time the average temperature of the loop 
increases (Figure |6]) while the maximum temper- 
ature (with thermal conduction included) stays 
approximately the same (Figure [5]). The reason 
for this is that as the current sheets evolve, mul- 
tiple reconnection acts as a continuous heating 
source to the plasma inside the loop, with the 
location of the heating changin g as the current 
sheets change ([Hood et al. I l2009h . 

The temperature response function for the 
SDO/AIA 171 A band is presented in Figure 
El The emissions in 171 A for SDO/AIA and 



TRACE are essentially the same, because their 
temperature response functions lie so close to- 
gether. The only difference between the two ob- 
servational platforms is the time resolution. With 
the TRACE exposure time 11 times longer than 
that of SDO/AIA, the images from SDO/AIA are 
sharper. 

In addition to the figures presented in this 
paper, the response functions from TRACE 195 
A and 284 A as well as those for the coronal band- 
widths of SDO/AIA were used in calculating the 
line of sight integral In all cases a broad 

range of temperatures are sampled from the sim- 
ulated data, with a large part of the temperature 
range common to all the response functions. As 
a result, the differences between the line inten- 
sity plots were minor; the same basic internal loop 
structure was observed with no new information to 
be gained from them. 

4. Conclusion 

The evolution of a coronal loop is studied, solv- 
ing the resistive MHD equations for a fully ionised 
plasma and with parallel thermal conduction in- 
cluded in the model. The loop is initialised as a 
straight cylinder with a twist above the stability 
threshold, which leads to the kink instability. Line 
of sight emission intensities were calculated of the 
simulation data, using the temperature response 
functions from TRACE. 

The simulations were initialised with physical 
parameters extracted from the observations of a 
coron al loop show n in Figure [1] ( Srivastaya et all 
2010h . Figure 5 in ISrivastava et al.l (|2010l ) shows 
the time evolution of a helically twisted structure 
of a flaring loop in TRACE 171 A that is observed 
for 4 minutes. The initial magnetic field structure 
used in the simulations is a kink-unstable force- 
free equilibrium where the twist varies with ra- 
dius (|Hood et al.ll2009h . This field structure was 
used in previous studie s of the kink instability 
(see iBotha et al.l (| 2 1 lh and references therein) 
and although its maximum twist is similar to the 
global estimated twist from the observations by 
Srivastava et al. ( 2010h . there is no guarantee that 
it is similar to the magnetic field of the observed 
loop. Throughout the evolution of the kink insta- 
bility, the internal structure of the simulated loop 
is shown in the generated images (Figure [7]) . This 
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structure evolves into a simpler configuration as 
the kink instability causes multiple reconnection 
events, which have the effect of straightening the 
internal magnetic field. I n con trast, the observa- 
tions bv lSrivastava et all ( 201Clh show a field struc- 
ture that stays largely intact for the duration of 
the observations. 

Footpoint brightening due to compression of 
the plasma is observed in the simulation re- 
sults, but is absent from the observations by 



Srivastava et al 



(2010). In the simulation this 
is due to the impenetrable boundaries of the nu- 
merical domain, which cause the plasma density 
to increase at the top and bottom boundaries of 
the loop axis (Figure On the Sun plasma 

moves through the transition region to the chro- 
mosphere and photosphere. It may be that foot- 
point brightening occurs in these lower regions, 
but the plasma is so cool that most likely it will 
not be captured by the response functions for 
coronal temperatures in spite of collectiv e foot 



study, which should be included at a later stage. 

AKS thanks Shobhna Srivastava for patient en- 
couragements. 
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Fig. 1. — The observed loop system used as source 
for the initialisation of the simulations. Half the 
loop length is shown with the helicity of right 
handed twist clearly visible. A time evolution of 
th e complete loop system is presented by Figure 5 
in lSrivastava et al] ( 2010( ). The image is in reverse 
colour. 
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Fig. 2. — Cartoon showing the orientation of the 
cylindrical loop in the Cartesian numerical do- 
main. The dotted line indicates the integration 
path of the line of sight integral ^ along the y 
direction and perpendicular to the (x, z) plane. 
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Fig. 3. — Normalised temperature response func- 
tions for TRACE 171 A (solid line) and AIA 171 
A (broken line). The two vertical dotted lines are 
respectively the minimum and maximum temper- 
atures during the simulation with parallel thermal 
conduction included, as measured in Figures [5] and 
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Fig. 4. — Intensity from the line-of-sight integration ([3]), using the simulation data as filtered through the 
TRACE 171 A temperature response function in Figure |3] (a) the simulation data at 417.6 s; (b) the time 
averaged data set; (c) the time and spatially averaged data set. The time exposure is 31.9 s and the spatial 
resolution is 0.375 Mm per pixel. The images are in reverse colour, with white representing the lowest value 
on the scale. The minimum and the maximum values of the reverse colour scale are the same as those used 
in Figure [7J These data sets were taken from the simulation with parallel thermal conduction. 
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Fig. 5. — Maximum temperatures during the evo- 
lution of the kink instability, for the two numerical 
runs with thermal conductivity included (broken 
line) and without thermal conduction (solid line). 
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Fig. 6. — Temperatures during the evolution of the 
kink instability for the two numerical runs with 
and without thermal conductivity. The mean tem- 
perature is represented by a solid line for the sim- 
ulation without thermal conduction and a broken 
line with thermal conductivity, similar to FigureJS] 
The mean temperature is calculated by averaging 
over the entire data cube. The minimum temper- 
atures are also included, with the dash-dot-dashed 
line the simulation without thermal conduction 
and the dotted line with thermal conductivity. 
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Fig. 7. — Simulated intensities as would be seen through the TRACE 171 A filter of the kink and its 
aftermath, with thermal conduction included in the calculation. The first image is taken at 290.0 s and the 
cadence is 31.9 s. The time exposure, spatial resolution and light intensity scale are the same as in Figure 
|4j The reverse colour table is such that the minimum intensity (white) is chosen so that the lowest 10% of 
the simulation values at time 290 s are eliminated from these plots. All the exposures use the same reverse 
colour scale. 
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Fig. 8. — Simulated intensities as would be seen through the TRACE 171 A filter of the kink and its 
aftermath, without thermal conduction. The exposure time, spatial resolution and light intensity scale are 
the same as in Figure [4j with the minimum intensity (white) and the maximum (black) of the reverse colour 
table having the same values as in Figure [7j The images are taken at the same times as those in Figure [7j 
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Fig. 9. — Temperature and mass density profiles 
along the length of the loop on its central axis at 
time 579.98 s, obtained from the simulation with 
thermal conduction. The solid line is the tempera- 
ture and the dashed line the density. Initialisation 
is with a constant temperature of 0.125 MK and 
a constant density of 1.67 ng m -3 . 
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